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Optically detected magnetic resonance (ODMR) is a way to characterize the NV^ centers. Recently, a re¬ 
markably sharp dip was observed in the ODMR with a high-density ensemble of NV centers, and this was 
reproduced by a theoretical model in [Zhu et al.. Nature Communications 5, 3424 (2014)], showing that the dip 
is a consequence of the spin-1 properties of the NV~ centers. Here, we present much more details of analysis 
to show how this model can be applied to investigate the properties of the NV^ centers. By using our model, 
we have reproduced the ODMR with and without applied magnetic fields. Also, we theoretically investigate 
how the ODMR is affected by the typical parameters of the ensemble NV“ centers such as strain distributions, 
inhomogeneous magnetic fields, and homogeneous broadening width. Our model could provide a way to esti¬ 
mate these parameters from the ODMR, which would be crucial to realize diamond-based quantum information 
processing. 


I. INTRODUCTION 

A nitrogen-vacancy (NV“) center in diamond M is a 
promising candidate to realize quantum information process¬ 
ing iil] and network lEl. An NV center is known to have 
a long coherence time such as a second llTll - [l3|] . The oper¬ 
ations such as qubit gates and measurements, which are ba¬ 
sic tools for quantum applications, have been demonstrated 
with a single NV center lfl4|] . Also, the entanglement gener¬ 
ation between distant nodes, which plays an essential role of 
quantum repeater, has been demonstrated by using photons as 
flying qubits emitted through distant two single NV centers 
11^ . An NV“ center can be used for a sensitive magnetic 
field sensor lIMi. An ensemble of NV centers can be 
also used for demonstrating quantum metrology ifTgl - l^ and 
physi cal p henomena in fundamental phy sics, such as quantum 
walk and quantum simulation 11^ . Also, the ensemble 
of NV“ centers can be used as the hybrid devices between 
different phy sical systems, in particular, superconducting sys¬ 
tems ll24l - f^ . Due to the effect of a superradiance, the ensem¬ 
ble of NV“ centers has a much stronger magnetic coupling 
with other systems than a single NV“ center. 

An NV“ center consists of a nitrogen atom and a yacancy 
in the adjacent site id, and this is a spin-1 system with three 
states of |0), I — 1), and |1). With a strong external magnetic 
field, the two exited states |1) and | — 1) of the NV“ center is 
energetically separated far from each other. In this case, the 
NV“ center can be considered as a spin ^ system by using 
a frequency selectiyity where |0) and |1) (| — 1)) constitute 
a qubit. On the other hand, with zero or weak applied mag - 
netic field, the NV“ center reyeals spin-1 properties 
Optically detected magnetic resonance (ODMR) is the gen¬ 
eral technique to inyestigate the properties of the NV“ centers 
13] . After applying a microwaye pulse, the NV“ centers are 
measured by an optical detection. Resonance obseryed with 
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specific microwaye frequencies let us know an energy struc¬ 
ture of a ground-state manifold of the NV“ centers. Also, we 
can estimate coherence properties of the NV“ center from the 
width of the peaks. 

Recently, a remarkably sharp dip has been obseryed around 
2870 MHz in the ODMR with zero applied magnetic fields 
ill m Ei- Although the ODMR results are usually fit 
by a sum of Lorentzians, the ODMR results obseryed in 
cannot be well reproduced by such a fitting 12^ . 
and no theoretical model can explain the dip until a new ap¬ 
proach is suggested in S. The model described in S con¬ 
tains spin-1 properties of the NV“ centers while most of the 
preyious models assume the NV“ center to be a spin-half sys¬ 
tem or use just a sum of Lorentzians to include the effect of 
the spin-1 properties in a phenomenological way a . By in¬ 
cluding the strain distributions, randomized magnetic fields, 
and homogeneous width of the NV“ centers, the sharp dip 
in the ODMR has been reproduced in 1^ . This model pro- 
yides us with an efficient tool to characterize the high-density 
ensemble of NV“ centers, which would be crucial to realize 
diamond-based quantum information processing. Moreoyer, 
this dip is shown to be the cause of a long-liyed collectiye 
dark state obseryed in a spectroscopy of superconductor dia¬ 
mond hybrid system, and so this dip could be useful if we will 
use the collectiye dark state for a long-liyed quantum memory 
of a superconducting qubit ifsill . 

In this paper, we present the details about how the model 
suggested in can be applied to inyestigate the properties 
of an ensemble of NV“ centers. An ensemble of NV“ cen¬ 
ters is affected by inhomogeneous magnetic fields, inhomoge¬ 
neous strain distributions, and homogeneous broadening. By 
taking into account of these as parameters in our model, we 
haye reproduced the ODMR with and without applied mag¬ 
netic field. Moreoyer, from a numerical simulation, we haye 
inyestigated how these parameters affect the sharp dip around 
2870 MHZ and the width of the each peak in the ODMR. 
We haye found that homogeneous broadening is releyant to 
change the dip in the ODMR. Also, we haye confirmed that 
the width of the peaks in the ODMR is insensitiye against the 
strain yariations if an external magnetic field is applied. More- 
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over, we have shown how our model could be used to estimate 
these parameters of the NV“ centers from the ODMR. 

The rest of this paper is organized as follows. In section 2, 
we explain the experimental setup. In section 3, we introduce 
the theoretical model introduced in JH.In section 4, we show 
the ODMR results and explain how these experimental results 
can be reproduced by our theoretical model. Finally, section 
5 contains a summary of our results. 


II. EXPERIMENTAL SETUP 


We begin by describing how we generate the NV“ centers 
in diamond. To create the NV“ center ensemble, we per¬ 
formed ion implantation of and we annealed the sample 

in high vacuum dHt]. The density of the NV“ centers is ap¬ 
proximately 5 X 10^^ cm“^, and we have the NV“ centers 
over the depth of Ipm from the surface of the diamond. 
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FIG. 1: NV“ center consists of a nitrogen atom (N) and a vacancy 
(V) in the adjacent site. Since NV” center is a spin-1 system, we 
have three states of |0), |1), and | — 1). We can characterize the 
NV” center by an optically detected magnetic resonance (ODMR) 
spectrum, and we perform the ODMR with an applied magnetic field 
of B = 0,1, 2 mT, along the [111] direction. 

The ODMR was performed on the diamond sample by a 
confocal microsco pe w ith a magnetic resonance system at 
room temperature lUin . We manipulate pulsed optical laser 
(532nm) and microwave independently. The magnetic field of 
0, 1, or 2 mT was applied along the [111] axis. With zero or 
weak applied magnetic field, a quantization axis of the NV” 
center is determined by the direction from the vacancy to the 
nitrogen, which we call an NV” axis. This axis is along one 
of four possible crystallographic axes. The NV” centers usu¬ 
ally occupy these four directions equally. The applied mag¬ 
netic field along [111] is aligned with one of these four axes 
as shown in Fig. [T] In this case, the Zeeman splitting of the 
NV” centers having the NV” axis of [111] is larger than that 
of the NV” centers having the other three NV” axes. 


III. MODEL 

We describe the model to simulate the ODMR of the NV” 
center ensemble, which was introduced in iS-TheNV axis 
provides us with the z axis. Microwave pulses are applied on 


the NV” centers, and the microwave pulses orthogonal to the 
z axis induce the excitation of the NV” centers. We define the 
X axis as such a orthogonal direction of the applied microwave 
at each NV” center. The Hamiltonian of the NV” centers is 
as follows. 


N 

H = [DkSl^ + - Si,) 

k^l 
( 


+ [Sx^kSy^k + Sy^j^Sx^k) + Xcos{Lot)S^ 
+ A^^Sz^klz.k H-,/c “t“ 


where S, (Ik) denotes a spin-1 operator of fcth elec- 

(k) 

tron (nuclear) spin, Dk denotes a zero-field splitting, E) 

(k) (k) 

(E 2 ') denotes a strain along x(y) direction, ge^^sBz • 

(— Iz,k) denotes a Zeeman term of the fcth electron 
(nuclear) spin, A denotes a microwave amplitude, oj denotes 
a microwave frequency, P denotes the quadrupole splitting, 
and yl|| (A_l) denotes a parallel (orthogonal) hyperfine cou¬ 
pling. For simplicity, we assume a homogeneous microwave 
amplitude Afc cx A. (In the appendix, we relax this constraint.) 
It is worth mentioning that the x and y component of the mag¬ 
netic field is insignificant to change quantized axis and so we 
consider only the effect of z axis of the magnetic field. Since 
the energy of the nuclear spin is detuned from the energy of 
the electron spin, the flip-flop term ^{S+^kI-,k + S-^kI+,k) 
is negligible and the parallel term A\\Sz,kIz,k is dominant. In 
this case, the effect of the nuclear spin is considered as ran¬ 
domized magnetic fields on the electron spin ll^[3^ . and we 
use this approximation throughout the paper. 

In a rotating frame defined by [/ = pgj._ 

form the rotating wave approximation, and we obtain the sim¬ 
plified Hamiltonian. 


N 

H^hJ2{iDk- u;)Sl, + eI^\sI, - Si,) 

k^l 

+E^l'\s,,kSy,k + Sy,,S,,,) + getiBBi'^'fSz + 

If the number of excitations in the spin ensemble is much 
smaller than the number of spins, we can consider the spin en¬ 
semble as a number of harmonic oscillators. In this case, we 
can replace the spin ladder operators as creational operators of 
the harmonic oscillators such as b\ ~ |H)fe(0|, dj, ~ |iA)fc(0| 
where |B)fc = ^(|l)fc +1 - l)fc), = ;^(|l)fc -1 -1^-)- 
By using this approximation, we can simplify the Hamiltonian 
as follows S. 

N 

H ^ {(wf ^ - uj)blb, + - uj)dldk 

k^l 

+ Jk{b\dk + hdi) + iJk{b\d, - bkd\) + -^{bk + 
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where = Dk + E['^\ Jk = g^iBB‘i^\ 

and J'f, = . 

The inhomogeneous broadening can be included in this 
model as following. We use Lorentzian distributions to in- 
elude an inhomogeneous effect of E^ , and E 2 [k = 
1, 2, • • • ,N). It is worth mentioning that the Lorentzian dis¬ 
tributions have been typically used to describe the inhomo¬ 
geneous broadening of the NV“ centers ll33l - l^ IddTl . For an 
inhomogeneous magnetic field Bz , we need to consider the 
following two effect. First, since there is an electron spin- 
half bath in the environment due to the substitutional N (PI) 
centers, NV“ centers are affected by randomized magnetic 
fields. Second, a hyperfine coupling of the nitrogen nuclear 
spin splits the energy of the NV“ center into three levels. So 
we use a random distribution of the magnetic fields with the 
form of the mixture of three Lorentzian functions. Here, each 
peak of the Lorentzian is separated with 27r x 2.3 MHz that 
corresponds to the hyperfine interaction with nuclear spin 

11^ 1^ . It is worth mentioning that, since the frequency shift 
of Dk is almost two-orders of magnitude smaller than that of 
E^'^ and E^'^ 1^ . we consider the effect of inhomogeneity 
of Dk as this order in this paper. 

We can describe the dynamics of the NV“ centers by using 
the Heisenberg equation as follows. 


-iYf,)hk-iJkdk +J'kdk-iX 
-iY^)dk-iJkbk-J'kbk (1) 


where rb(= F^) denotes the homogeneous width of the NV 
center. We assume that the initial state is a vacuum state. 
Since we consider a steady state after a long time, we can 
set the time derivative as zero. In this condition, we obtain 

(^l,t=oo^fc.t=oo) 

^ I_ \{uj - + iVd) _,2^2) 

(u; - + tF,) - {\Jk? + 14 P) 

(dkd=oodk,t=oo) 

_ I_ A(Jfc — iJ'k) _,2^2^ 

(w - -I- irb)(a; - + iTd) - (|Tfcp + I4P) 


The average probability of the NV center to be in the energy 
eigenstates other than |0) can be calculated as 

1 ^ . 

+ i‘^k,t=oodk,t=co))- ( 4 ) 

k=l 


In the actual experiment, if we excite the NV“ centers by 
the microwave pulses, the intensity of the photons emitted 
from the NV“ centers will be changed from the baseline emis¬ 
sion rate Iq. This change is linear with Pe- So, to fit the ex¬ 
periment with our model, we use a function of {Iq — aP^)lI q 
where a denotes a fitting parameter , and this corresponds to 
the ODMR signals. 


IV. MAIN RESULTS 
A. Reproducing the experimental results 



FIG. 2: ODMR with zero applied magnetic fields. S{gfj,BBz)/2n = 
1.96 MHz (HWHM), SEi/2-k = SE2/2tt = 0.73 MHz (HWHM), 
SDk = 0.01 MHz (HWHM), X/2tv = 2 MHz, r6/27r = rd/27r = 
0.3 MHz. Also, we assume a Nitrogen hyperfine coupling of 27r x 2.3 
MHz. The red line denotes a numerical simulation and blue dots 
denote the experimental results. 



Frequency (MHz) 


FIG. 3: ODMR with an applied magnetic field of 1 mT. We use 
the same parameters as the Fig. The red line denotes a numerical 
simulation and blue dots denote the experimental results. 

By using the model described above, we have reproduced 
the ODMR signals when we apply B = 0,1, 2 mT, as shown 
in Figs |2] 0 and |4] A sharp dip is observed around 2870 
MHz for the case of i? = 0 mT, and our simulation can re¬ 
produce this. Two peaks are observed in the ODMR with 
zero applied magnetic field as shown in Fig |2] which corre¬ 
sponds to the transition between the state |0) and one of the 
other energy eigenstates. If we consider a single NV“ cen¬ 
ter, the frequency difference between the two exited states is 
Suj = 2yJ{ggBBzy + {EiY -F Since we consider an 

ensemble of the NV“ center, this frequency difference vari¬ 
ates depending on the position of the NV“ center. For sim¬ 
plicity, we use a dimensionless variable for Bz, Ei, and E 2 , 
defined as Bz = gfJ-BBz/^, Ei = Eij^, and E 2 = ^^ 2/7 
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Frequency (MFIz) 

FIG. 4: ODMR with an applied magnetic field of 2 mT. We use 
the same parameters as the Fig. The red line denotes a numerical 
simulation and blue dots denote the experimental results. 

where 7 denotes a damping rate with an unit of the frequency. 
To calculate a probability that the two energy eigenstates such 
as |1) and | — 1) are degenerate (Suj = 0), we define probabil¬ 
ity density functions of Ei, and E 2 as Pa[B^), Ph{Ei), 
and Pc{E 2 ), respectively. The joint probability is calculated 
as 

P{Ezt El, i?2)Ai?^A.Ei A.E 2 
= Pa{B, = 0)AB, ■ PkiEi = 0)AEi ■ P,{E2 = 0)AE2. 

where AB^, AEi, and AE 2 denote a finite range of each vari¬ 
able. We assume = Pa{B,)Pb{Ei)P,{E 2 ) 

because these are independent. By using spherical coordinates 
where B^ = r sin 0 cos (j), Ei = r sin 9 sin cj), E 2 = r cos 9 

with r = \J\Bz\'^ + \Ei\'^ + we rewrite this as 

P{Bz;, El, E2 )ABzAEiAE 2 
= P^{B^ = 0)PbiEi = 0 )Pc{E2 = sin6»ArA6lA<^. 

This shows that, even if we consider a finite range ABz, AEi, 
and AE 2 , the probability for the two energy eigenstates to be 
exactly degenerate (r = 0) is zero. This means that, if homo¬ 
geneous broadening is negligible, the two peaks to denote the 
two energy eigenstates of each NV“ center should be always 
separated in the ODMR so that the ODMR signal at the fre¬ 
quency of £)/27r = 2870 MHz should be the same as the base 
line. However, due to the effect of the homogeneous broaden¬ 
ing, small signals deviated from the base line can be observed 
at the frequency of I?/27r = 2870 MHz. This is the cause 
of the sharp dip observed around the frequency of 27r x 2870 
MHz in the ODMR. 

With an applied magnetic field, four peaks are observed in 
the ODMR where two of them are larger than the other two, 
as shown in Figs |3] and |4] The two smaller peaks correspond 
to the energy eigenstates of the NV“ centers with an NV“ axe 
along [111], which is aligned with the applied magnetic field. 
A quarter of the NV“ centers in the ensemble have such an 
NV“ axis. The other larger peaks come from the other NV“ 
centers where the applied magnetic field is not aligned with 


the NV“ axis. Three-quarters NV“ center have such axes. 
In this case, the Zeeman splitting of these is smaller than that 
of the NV“ centers with the [111] axis. It is worth mention¬ 
ing that a small dip is observed in the ImT ODMR around 
27 r X 2870 MHz due to the mechanism explained above. On 
the other hand, such a dip is not clearly observed in the 2mT 
ODMR, because the NV“ centers are considered to be as ap¬ 
proximate two-level systems in this regime. 


B. The behavior of the ODMR against the change in the 
parameters 

We perform a numerical simulation with several parameters 
to understand the behavior of the sharp dip. In the Fig|5]a, we 
change the parameter Ff, while we fix the other parameters. 
Similarly, in the Fig |5] b (c), we change the parameter 5Bk 
(SEk) while we fix the other parameters. We have found that 
the sharp dip is very sensitive against the change in Fwhile 
the dip is relatively insensitive against the change in SEk and 
SBk. 

Also, we perform a numerical simulation with several pa¬ 
rameters for the ODMR with an applied magnetic field. In the 
Figs |6l we have plotted one of the peaks of the ODMR with 
an applied magnetic field of 2mT. This peak corresponds to 
a transition between |0) and | — 1) of the NV“ center with 
an axis of [111]. From the numerical simulations, we have 
found that this ODMR signals with applied magnetic field is 
robust against the strain variations SEk, while the peak will 
be broadened due to the effect of the randomized magnetic 
field SBk- The frequency difference between the ground state 
and another energy eigenstate can be calculated as Soj' = 

Dk — \J E^'^ + E^'^ + If the applied magnetic 

field is large, we obtain Soj' ~ ± ^ ^ ■ 

This means that the effect of the strain is insignificant in this 
regime while the inhomogeneous magnetic field from the en¬ 
vironment can easily change this frequency. These can ex¬ 
plain the simulation results shown in Figs |6] where the change 
of inhomogeneous magnetic fields affects the width of the 
peak while the peak width is insensitive against the inhomoge¬ 
neous strain. Such an effect to suppress the strain distributions 
by an applied magnetic field was mentioned in and was 
recently demonstrated in a vacuum Rabi oscillation between 
a superconducting flux qubit and NV“ centers in 1^ . Our 
results here are consistent with these previous results. 


C. Parameter estimation 

An ensemble of NV“ centers is affected by inhomogeneous 
magnetic fields, inhomogeneous strain distributions, and ho¬ 
mogeneous broadening. In the ODMR, the observed peaks 
contain the information of the total width that is a compos¬ 
ite effect of three noise mentioned above, and so it was not 
straightforward to separate these three effects for the estima¬ 
tion about how individual noise contributes to the width. 
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FIG. 6: Numerical simulation of ODMR with an applied magnetic 
field of 2 mT. Here, x axis denotes the microwave frequency, y axis 
denotes SEl2n for the left figure while y axis denotes S^gysB)/2'n 
for the right figure, and 2 : axis denotes the ODMR signal intensity. 
Other than SE or S{giJ.BBz)/2n, we use the same parameters as 
the Fig. [4] These show that the peak width is much more sensi¬ 
tive against the inhomogeneous magnetic fields than inhomogeneous 
strain. 


is much smaller than the 6Bk and 6Ek ll^ 1^ . and 

so it seems that the effect of Fb might be hindered by a huge 
influence of SBk and SE^- However, since the dip is sensitive 
against the change in Fb, we could accurately estimate the Fb 
even under the effect of 5Bk and SE^- Secondly, as we have 
shown, the ODMR signals with applied magnetic field is ro¬ 
bust against the strain variations 6Ek while the peak will be 
broadened due to the effect of the randomized magnetic field 

(k) 

SBk- We can use these properties to estimate the 6Bl ' and 
SEk- Since the ODMR with an applied magnetic field is in¬ 
sensitive against SEk, we can estimate SBi^^ from this exper- 
imental data. Since we have estimated SBz ' and Fb from the 
prescription described above, we fix these parameters so that 
we can estimate SEk from the ODMR with zero applied mag¬ 
netic field. Therefore, by applying these procedure, we could 
estimate the parameters of the NV“ centers such as inhomo¬ 
geneous magnetic fields, inhomogeneous strain distributions, 
and homogeneous broadening. 


V. SUMMARY 


FIG. 5: Numerical simulation of ODMR with zero applied magnetic 
fields. Here, x axis denotes the microwave frequency, y axis denotes 
7/27r (for the figure a) or <5 {gyBBz)/2n (for the figure b) or SE /27r 
(for the figure c), and z axis denotes the ODMR signal intensity. 
Other than the inhomogeneous width of the inhomogeneous width, 
we use the same parameters as the Fig. 


Interestingly, our model could be used to determine these 
three parameters by reproducing the ODMR with and without 
applied magnetic fields. Firstly, as we described before, the 
sharp dip in the ODMR is very sensitive against the change in 
Fb, while the dip is relatively insensitive against the change 
in SEk and SBk- These properties are important to determine 
the value of Fb from the analysis of the ODMR. Usually, Fb 


In conclusion, we have studied an ODMR with a high- 
density ensemble of NV“ centers. Our model succeeds to re¬ 
produce the ODMR with and without applied magnetic field. 
Also, we have shown that our model is useful to determine 
the typical parameters of the ensemble NV“ centers such as 
strain distributions, inhomogeneous magnetic fields, and ho¬ 
mogeneous broadening width. Such a parameter estimation is 
essential for the use of NV“ centers to realize diamond-based 
quantum information processing. 

Y.M thanks K. Nemoto and H. Nakano for discussion. This 
work was supported by JSPS KAKENHI No. 15K17732, 
JSPS KAKENHI Grant No. 25220601, and the Commis¬ 
sioned Research of NICT. 
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VI. APPENDIX 

Here, we consider the effect of inhomogeneous microwave 
amplitude. If we have such an inhomogeneity, by solving the 
Heisenberg equation, we obtain 

^\,t=c<fik,t=oo) 

I_ Xkjuj-ujf'' + trd) _2 

(c. - ccf + *r,)(a; - tuf + *r,) - (| + | J'P) 

{^k,t=oo^k,t=oo) 

_ I_ XkjJk — iJ'k) _12 

" ' (u; - + *r,)(cc - + zFd) - (I JfcP + 14P) 

where the value of A differs depending on the position of the 
NV“ centers. If we define an average probability of the NV“ 
center in the bright (dark) state as Pi, {Pd), we obtain 


m flooK^) 

Yfi floor(-^) 

= E of) m 

j^l ^m' k^l 

where 

(fe) ^ I_ {uj-J,i^^+iTd) _ 2 

^ (w - wf’ + iTb){uj - Jp + ikd) - (I JfcP + I4P) 

(d) ^ I_ (Jfc - iJ'i,) _,2 

^ (w +irb){uj-u}P +ird) - (|Jfep + I4P) 

Therefore, we obtain 


N 




I__I 


k=l 


Pd = 




N 


I _ >^k(Jk-iJi) _12 

I (^^p)+irb)(<..;-<..;‘'“4ird)-(|JfcP + | J(P) ' 


N 


(5) 


Since inhomogeneity of A is independent from the inhomo¬ 
geneity of LtJb, u!d, J, and J', we can rewrite these probabilities 
for a large number of NV“ centers as follows 


N' 

k^l 

1 

Pd ^ ^9) 

where |Aavp = and iV' = floor(^). The 

probability of the NV“ center in the ground states can be cal¬ 
culated as 


Pe — Pb + Pd 


( 10 ) 


Pb 



i=i 


floor(i^) 

E 


fc=i 



and this is the same form as the probability of the homoge¬ 
neous microwave amplitude case described in the Eq. (|4| 
where N (A^) is replaced by N' (jAp^). So the inhomoge¬ 
neous microwave amplitude does not affect the theoretical 
(6) prediction of ODMR signals. 
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